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ABSTRACT 



We have selected a small sample of post-AGB stars in transition towards the planetary 
nebula and present new Very Large Array multi-frequency high-angular resolution radio ob- 
servations of them. The multi-frequency data are used to create and model the targets' radio 
continuum spectra, proving that these stars started their evolution as very young planetary neb- 
ulae. In the optically thin range, the slopes are compatible with the expected spectral index 
(-0.1). Two targets (IRAS 18062+2410 and 17423-1755) seem to be optically thick even at 
high frequency, as observed in a handful of other post-AGB stars in the literature, while a third 
one (IRAS 20462+3416) shows a possible contribution from cold dust. In IRAS 18062+2410, 
where we have three observations spanning a period of four years, we detect an increase in its 
flux density, similar to that observed in CRL 618. 

High-angular resolution imaging shows bipolar structures that may be due to circumstel- 
lar tori, although a different hypothesis (i.e., jets) could also explain the observations. Further 
observations and monitoring of these sources will enable us to test the current evolutionary 
models of planetary nebulae. 
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1 INTRODUCTION 

The evolution of planetary nebulae remains a challenging topic in 
astrophysics. Several studies have focused on the link between the 
post-Asymptotic Giant Branch (AGB) and planetary nebula phases 
to investigate the changes that occur in the evolutionary stage be- 
tween these object classes. Planetary nebulae (PN) have been ob- 
served over a wide wavelength range, from X-ray to radio frequen- 
cies. Their complex morphologies and the shaping mechanisms that 
produced them are still a matter of debate. Companion stars, jets 
from central stars, magnetic fields, dust tori, and interacting winds 
are some of the possible shaping agents suggested as being respon- 
sible for various PN morphologies, and an overlap of their actions 
cannot be ruled out. 

Hubble Space Telescope (HST) imaging has demonstrated 
that pre-PN show non-spherical symmetry even at late spectral 
types, which implies that th e shaping process in t hese sources starts 
very early after the AGB (Balick & Frank 2002). Non-spherically 
symmetric structures hav e also been observed i n AGB stars that are 
part of a binary system jKarovska et~aT]|2007h . which could indi- 
cate that some shaping agents are already active in this evolution- 
ary phase. In general, the current theory of PN evolution is based on 
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the Interacting Stellar Wind (ISW) model iKwok et alj| 19781) and 
its generalised version (Kah n & Westlll985h . However, this model 
does not account for shaping because it assumes that an asymmetric 
distribution of matter is already present when the wind interaction 
occurs. Other models take into account the possibl e role of jets, and 
have been successfully applied to some nebulae ( Sah ai & Trauger] 
1998). It has been suggested that the interaction with a companion 
obj ect, even a m assive planet, may provide the necessary asymme- 
try (Soker 2006). Also, large scale magnetic fields, influencing or 
determining the shape s, might be sustained by a dynamo process 
jBlackman et all200lh . 

The birth of a PN is defined by the observation of an ioni- 
sation front, which is itself a shaping agent and could heavily in- 
fluence the morphology established in earlier evolutionary phases, 
disrupting the molecular and dust circumstellar shells. Near- and 
mid-IR images of PN, which trace the molecular and warm dust 
emission, have been compared to optical line images, which trace 
the io nised elements, an d have shown the presence of similar struc- 
tures dLatteret all 1 19951) . Therefore it seems that molecular and 
ionised gas and dust grains can spatially coexist in these sources. 
This makes the investigation of the spatial distribution and physical 
properties of the ionised component in these envelopes even more 
compelling. 

In this context, important information can be provided by ob- 
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servations of very young PN, or pre-PN, where the physical pro- 
cesses associated with PN formation are still occurring. To investi- 
gate the properties of these rare objects in transition from the post- 
AGB to PN, we have selected a sample of pre-PN and searched for 
radio emission from ionised shells. The targets were selected from 
stars classified in the literature as hot post-AGB candidates, show- 
ing strong far-IR excess and B spectral type features. We detected 
radio emission in 10 sources in the selected sample dUmana et al.l 
120041) . The detection of radio continuum emission is proof of the 
presence of free electrons and therefore of an ionised shell. 

For a better understanding of these sources, we have per- 
formed a multi-frequency follow-up to build up radio continuum 
spectra and collect information about the physical conditions in 
the nebulae. Several nebular parameters can be determined by ra- 
dio flux density measurements (for example, electron density and 
ionised mass) but most calculations rely on the assumption of an 
optically thin emission. Although this is presumably the case for 
observations at frequencies higher than 5 GHz, only the construc- 
tion of the whole centimetre continuum spectrum can enable us to 
distinguish between optically thick and thin regimes. VLA A array 
observations have also given us insight in the morphologies of the 
envelopes. 



2 OBSERVATIONS AND RESULTS 

2.1 Low-angular resolution observations 

Multi-frequency observations were carried out at the Very Large 
Array, managed by NRAO, in several runs, at 1.4, 4.8, 8.4, 14.9, 
and 22.4 GHz (respectively 20, 6, 3.6, 2, and 1.3 cm), in D and 
C arrays, to inspect the spectral energy distribution in the radio 
range. Table Q] lists the observati on details. The 8.4 GHz obser- 
vations were performed in 2001 dUmana et al]|2004l) . During the 
2003 runs, strong interference led to the failure of almost all of 
the 1.4 GHz observations and part of the 14.9 GHz ones. To ac- 
quire information about the lower frequency range of the spectrum 
of four sources (namely IRAS 18062+2410, 19336-0400, 19590- 
1249, 20462+3114), observations were performed in 2005 at both 
1.4 and 8.4 GHz; the latter ones enabled us to check for any in- 
trinsic variation in the emitted flux. Although these latter runs were 
carried out in a different array configuration (with a smaller synthe- 
sised beam), this cannot cause missing flux, since all of our sources 
are point-like in both D and C arrays at any frequency. The Largest 
Angular Scale (LASflis much larger than the typical sizes of these 
objects (^ 2") even at the shortest wavelength, being 60" at 22.4 
GHz (D and C arrays have the same LAS). The nominal beam sizes 
in D array are 44", 14", 8.4", 3.9", and 2.8", while in C array they 
are 12.5", 3.9", 2.3", 1.2", and 0.9", respectively at 1.4, 4.8, 8.4, 
14.9, and 22.4 GHz. One of the nebulae detected in 2001 (18371- 
3 159) could not be observed because of limited observation time. 

A typical observing scan consisted of two 15-minute point- 
ings on target, preceded and followed by 1 -minute scans on the 
phase calibrator. For a better phase calibration, 22.4 GHz observa- 
tions were carried out pointing at the target for four 7-minute scans 
alternating with 1 -minute scans on the phase calibrator. 

The data were reduced with the Astronomical Image 
Processing System (AIPS), according to the recommended reduc- 
tion process: the data set was FILLMed into AIPS with average 

1 The LAS is the size of the largest structure that an interferometer can 
map. 



opacity correction and nominal sensitivities (doweight=\), edited 
to flag interference or any other bad points and CALIBrated along 
with point weights (docalib=2). Maps were obtained using the task 
IMAGR with natural weights and CLEANed down to the theoreti- 
cal rms noise whenever possible, performing some hundreds of it- 
erations of the CLEAN algorithm. The flux density for each source 
in every map was estimated by fitting a Gaussian to the unresolved 
source (task JMFIT), and the rms noise was calculated in an area 
much larger than the synthesised beam (>100 beam), without evi- 
dent sources in it (task IMEAN). 

For IRAS 17423-1755 and 06556+1623, a second nearby 
(within 40") source was detected and large wavelength data cannot 
be considered indicative of their flux densities because of confusion 
problems. Tables Q] and [2] summarize the observation runs and the 
flux density measurements obtained. The reported error has been 
calculated as a — -J rms 2 + ((J ca iSi,) 2 , where a ca i is a 5% abso- 
lute calibration error. 



2.2 High-angular resolution observations 

To inspect the morphology of our targets, we carried out obser- 
vations at 8.4 GHz in A array at the VLA for those sources that 
were brighter tha n 1 mjy at this frequ ency in our detection exper- 
iment reported in lUmana etail d2004l) : IRAS 17381-1616, IRAS 
18062+2410, IRAS 18442-1144, IRAS 19336-0400, and IRAS 
19590-1249. The 1 mjy limit was set to select nebulae with sur- 
face brightnesses high enough to be imaged with the A array. The 
angular resolution of these observations was ~0.2". Two runs were 
performed in 2003: the first on July 10 (UT 07:04:50; 10:29:30) and 
the second one on August 15 (UT 05:12:30; 06:37:50). The inte- 
gration time per object was about 1.5 hours. The reduction process 
followed the standard format, as reported in the previous section. 
In the task IMAGR, the pixel size was set to 0.05", so that the con- 
volution beam is at least three pixels wide; Brigg's weighting of the 
data was used to find a good compromise between sensitivity and 
resolution (robust= 0) for all the resolved sources except 19590- 
1249, as natural weights were necessary in order not to miss its faint 
emission. For 17381-1616 and 18062+2410, we set respectively ro- 
bust= -2 and robust= -4 (nearly pure uniform weighting) to increase 
the resolution. The number of CLEAN iterations was set to 1000, 
minpatch 127, and gain 0.1. The absolute calibrator was 3C48, its 
flux density set to 3.15 Jy. On October 29 2004 we also managed 
to perform a snapshot (15 minutes on source) in A array of IRAS 
22023+5249, which we had detected as a radio source in a previ- 
ous observation run in C array and had measured as 2.0±0. 1 mjy 
(including absolute flux density calibration error). 

All the selected targets showed clear bipolar structures ex- 
cept IRAS 17381-1616 and IRAS 18062+2410, which were barely 
(IRAS 17381-1616) or not (IRAS 18062+2410) resolved. FigureQ] 
shows the structures found in all the radio envelopes as observed 
with the A array. 



3 RADIO PROPERTIES 
3.1 Modelling 

We have modelled the radio continuum spectra assuming that the 
central star is surrounded by a shell of ionised gas with density 
radially decreasing as r~ 2 . We have calculated a 3D distribution 
of mass in a spherical shell with outer radius Rout, inner radius 
Ki n = rjR ou t (with < r\ < 1), and density at the inner radius 
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Figure 1. Radio maps at 8.4 GHz of the targets observed in high angular resolution. The flux density levels, in units of /xjy/beam, are: in 17381 (-3, 3, 4, 5, 6, 
7, 8, 9)<r; in 18062 (-3, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30)cr; in 18442 (-3, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36)cr; in 19336 (-3, 3, 5, 7, 9, 11, 13, 15, 17, 
19, 21)<t; in 19590 (-3, 3, 4, 5, 6, 7, 8, 9, 10)<r; in 22023 (-3, 3, 4, 5, 7, 9, 1 1, 13)<r; where a indicates the rms calculated in each map and listed in Tabled] The 
maps do not show any negative level, because no such levels are present in the mapped areas. 
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Table 1. Summary of all the performed observation runs. 



Low-angular resolution 


Date 


VLA Conf. 


Freq. (GHz) Flux Calibr. 


Targets 


2003 
Feb 15 


D 


1.4,4.8, 8.4, 14.9, 22.4 


3C286 


18062+2410 
17423-1755 
17460-3114 
20462+3416 


Feb 16 


D 


1.4,4.8, 8.4, 14.9, 22.4 


3C286 


20462+3416 
19336-0400 
18442-1144 
19590-1249 


Mar 8 


D 


1 A A Q Q A 1/1Q A 
i.'r, 4.5, o.4, V s *. y, ZZ.4 




17381-1616 
17460-3114 
18062+2410 


Mar 21 


D 


1.4,4.8, 8.4, 14.9, 22.4 


3C48 


06556-1623 


2004 
Apr 18 


C 


8.4 


3C48 


22023+5249 


2005 
Jul 25 


C 


1.4 


3C48 


18062+2410 
19336-0400 
19590-1249 


Jul 26 


C 


1.4, 8.4 


3C48 


18062+2410 
19336-0400 
19590-1249 
20462+3416 


High-angular resolution 


2003 
Jul 10 


A 


8.4 


3C48 


18062+2410 
18442-1144 
19590-1249 
19336-0400 


Aug 15 


A 


8.4 


3C48 


17381-1616 


2004 
Oct 29 


A 


8.4 


3C286 


22023+5249 



pin- The calculation of the average density in the shell links this 
parameter to the density at the inner radius: < p >~ rj 2 (l~ri)p in . 
The electron temperature was set to 10 K for all nebu lae. Details 
of the radio model can be found in lUmanaet ai]( l2008h . 



Only a few estimates of the distances to the stars in our sam- 
ple are available. Where available in the literature, we have used the 
most common estimates of this parameter. For IRAS 18442-1244 
and 19336-0400, since their sizes and flux densities are comparable 
to what detected in IRAS 19590-1249, we assume the same dis- 
tance as for this source. For the remaining targets we have adopted 
a standard value of 1 kpc as an estimate of the order of magni- 
tude of their distances. Our radio model code requires as an input 
the distance to the star as well as the density at the inner radius of 
the shell, outer radius of the ionised shell, inner to outer radius ra- 
tio, and electron temperature. For the outer radii we have a ssumed 
as upper limits the sizes estimated in lUmana et al] d2004h . which 
were calculated from the VLA theoretical angular resolution. For 
those objects that were observed with the VLA in A array (IRAS 
17381-1616, 18062+2410, 18442-1144, 19336-0400, and 19590- 



12490, we used our estimate of their radii, based on their 3cr size. 
For IRAS 18062+2410 it was necessary to use a different geome- 
try, since with a shell structure we could not model a steep enough 
increase of flux density in the optically thick range. Therefore we 
used a sphere instead of a spherical shell. 

The results of our radio models are reported in Table [3] The 
density and ionised mass values are in agreement with those ex- 
pected for young planetary nebulae. For example, in NGC 7027, a 
high excitation y oung PN, values of 10 4 cm~ 3 and 0.005 Mq have 
been measured dBains et alj2003h , although, even in this very well 
studied PN, these results suffer from the distance uncertainty. Such 
distance independent parameters as emission measure and bright- 
ness temperature ha ve been previously calculated for our targets 
dUmana et al.l 120041) and found to match the expected values for 
young PN. 

We notice that in some cases the data points, especially at 



2 We do not have multi-frequency observations of IRAS 22023+5249, 
since we had not radio detected this target at the time of our multi-frequency 
run. 
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Table 2. Flux densities and their errors (including a 5% absolute flux density calibration error) from our D and C array observation runs are shown for each 
source and frequency. The asterisks indicate the sources chosen for high-angular resolution imaging. The spectral indexes in the optically thin range and the 
rms in the high angular resolution maps are also shown. 



IRAS ID 




Flux density (mJy) 






Spectral index 


A array rms (mJy) 


A array beam size (") 




1.4 Ghz 


4.8 Ghz 


8.4 GHz 


14.9 GHz 


22.4 GHz 




8.4 GHz 


8.4 GHz 


06556-1623 


2.0±0.1 1 


0.60±0.04 


0.55±0.03 


0.43±0.09 


0.36±0.06 


-0.29±0.20 






17381-1616* 


1.1±0.1 


1.61 ±0.09 


1.42±0.09 


1.3±0.1 


1.19±0.08 


-0.20±0.11 


0.038 


0.26 


17423-1755 


_2 


1.74±0.06 1 


0.26±0.03 


_2 


0.75±0.08 


+1.08 5 






17460-3114* 


_2 


1.6±0.1 


1.29±0.08 


1.7±0.1 


1.24±0.08 


-0.10±0.12 






18062+2410* 


0.26±0.04 3 


1.27±0.07 


1.46±0.09 4 


_2 


1.8±0.1 


+0.22±0.10 


0.045 


0.22 


18442-1144* 


_2 


17.8±0.9 


19.2±0.9 


14.2±0.7 


13.5±0.7 


-0.22±0.09 


0.038 


0.28 


19336-0400* 


8.0±0.4 3 


10.4 ±0.5 


9.7±0.5 


8.0±0.4 


7.0±0.3 


-0.26±0.09 


0.034 


0.27 


19590-1249* 


3.2±0.2 3 


3.1±0.2 


2.8±0.1 


2.3±0.1 


2.1 ±0.1 


-0.26±0.10 


0.028 


0.32 


20462+3416 


0.48±0.06 3 


0.60±0.05 


0.42±0.05 


0.50±0.09 


0.83±0.07 


-0.02±0.17 6 






22023+5249 






2.0±0.1 








0.027 


0.29 



1 More than one source in the beam: the flux density must be regarded as an upper limit. 

2 The data were not recovered. 

3 The measurement was performed in 2005. 

4 In 2005 the flux density was 2.1 ±0.1 mJy. 

5 This value is calculated with two data points only. 

6 The 22.4 GHz point was not included in the fitting procedure. 

Table 3. Model parameters for the radio shells. For 18062 a solid sphere geometrical model was used. IRAS 18062+2410 was modelled as an emitting sphere, 
not a shell. 



Target 


Distance 
kpc 


Rout 
II 


R-in /Rout 


Pin 

10 4 cm~ 3 


<P> 
10 3 cm" 3 


M ion 
10~ 4 M 


06556-1623 


1 


0.5 


0.15 


15 


3 


0.1 


17381-1616 


1 


0.35 


0.195 


30 


9 


0.14 


17460-3114 


1 


1.1 


0.049 


40 


0.9 


0.4 


18062+2410 


6.4 1 


0.061 




25 


8 


1.7 


18442-1144 


4 


1 


0.225 


8.8 


3.5 


77 


19336-0400 


4 


1 


0.2 


8 


2.6 


57 


19590-1249 


4 2 


1.2 


0.28 


2 


1.1 


40 


20462+3416 


3 3 


1.1 


0.11 


4 


0.04 


5 



^ BogdanovlfeqOOl) 

2 ConlonetalJJl993l) 

3 lArkhipova et alj fepoil) 

higher frequency, seem to show a steeper than expected slope in 
the optically thin range (Si, ~ u~ 01 ). Besides our modelling, we 
estimated the empirical spectral indexes by a linear fitting (in the 
log-log plane) of the data points that correspond to the optically 
thin regime (typically 4.8-22.4 GHz data). The results of the fitting 
are summarized in Table [2] 

Because at high frequency (14.9 and 22.4 GHz) a lack of short 
baselines may cause a lower sensitivity to flux from extended struc- 
tures (therefore a lower total flux density), we checked whether ta- 
pering our data (i.e., assigning larger weights to short baselines) 
could lead to larger flux densities, thus making the slopes in the 
spectra less steep. We did not find any significant difference in the 
results from tapered and non-tapered data. The empirical spectral 
indexes of three of the targets (IRAS 06556-1623, 17381-1616, and 
17460-3 1 14) are close to the expected value of -0.1, and out of the 
expected range for three others (IRAS 18442-1144, 19336-0400, 
and 19590-1249). Although we cannot rule out that the observed 
steeper slopes might be real and further observations should be car- 
ried out, if we consider that the reported errors are based on a statis- 
tics on 4 points only, the deviation from -0.1 is not statistically sig- 



nificant. Therefore we conclude that, within the errors, the spectral 
indexes are compatible with their expected value. 

Figure [2] shows the data points obtained in our VLA runs, 
the model spectra, and the linear fitting. We have not included 
IRAS 17423-1755 in our modelling, because only two data points 
are available. The flux density values obtained for this target would 
indicate an optically thick spectrum up to 22 GHz. 

Another particular result is the spectrum of 
IRAS 20462+3416, which is well matched by the model in 
the 1-15 GHz range, while the 22.4 GHz point is higher than the 
others. Since the measured radio flux densities are all below 1 mJy, 
even a small contribution from cold dust may be not negligible 
and so explain the 22.4 GHz measurement. Observations in the 
mm and sub-mm ranges would assess whether such a cold dust 
component is present or not. 



3.2 IRAS 18062+2410 

In our modelling we have used data from 2001, 2003, and 2005, 
since the sources that were re-observed in 2005 (see above) did not 



6 L. Cerrigone et al. 




Figure 2. Plots of the multi-frequency data obtained at the VLA. For 18062+2410, 19336-0400, 19590-1249, and 20462+3416, 2005 1.4 GHz data points have 
also been used. In 18062 they have been over-plotted as diamonds. The dashed lines are the result of the linear (in the log-log plane) fitting of the optically 
thin range points. The plots therefore show the difference between the modelled (solid line) and empirical (dashed line) dependence on frequency. The error 
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IRAS 18062+2410 

3.0 r — 1 1 1 1 ' ' 1 1 
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1.0 - 

o.5 r 

2000 2002 2004 2006 
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Figure 3. Increase of the 8.4 GHz flux density in IRAS 18062+2410 during 
years 2001-2005. 



be set in the standard ISW scenario, require a strong density gradi- 
ent in the nebula, with higher density regions p erpendicular to the 
main axis of the nebula JAaquist & Kwokl 19961) . Following this in- 
terpretation, the two bright peaks of emission are due to an opacity 
effect, with the optically thin radio flux tracing higher density re- 
gions. The other possibility is that we are actually observing the 
very early onset of outflows/jets piercing the ci rcumstellar enve- 
lope, as evident in many HST images of pre-PN (Sa hai & Traugej 
1998) and also pointed out by high-resoluti on CO observations o f 
the molecular component of the envelope (Hug gins etail feoOO). 
However, it is difficult to interpret the morphology observed in 
IRAS 22023+5349 in the framework of the standard ISW model, 
even invoking a strong density gradient in the nebula. For this ob- 
ject, a jet would be a more likely source of the observed morphol- 
ogy- 



show any flux density variation with respect to the 2001 8.4 GHz 
results within observational errors. The only exception is IRAS 
18062+2410, whose flux density varied by about +44%, increasing 
from 1.46±0.09 in 2001 to 2.1 ±0.1 mjy in 2005, while its sec- 
ondary calibrator, 1753+288, varied from 0.82±0.03 Jy in 2001 to 
1.41±0.04 Jy in 2005 (+75%). Unfortunately it was not possible to 
use the same primary calibrator in the two runs (see Table[TJ, which 
introduces a different density scale for the two measurements, al- 
though the error should be within a few percent (3-5%). The flux 
density variation in the secondary calibrator might indicate issues 
with the flux density calibration. Since other targets were also ob- 
served in the same run and no change was measured in their flux 
densities within the observational uncertainties, we conclude that 
the variation observed in 1753+288 is real, as possible for a sec- 
ondary calibrator, but it does not account for the variation observed 
in the target object. When we combine all of our observations of 
this target, including our 2003 measurement (A array data-set) of 
1.66±0.09 mJy, we can conclude that the radio flux density of 
IRAS 18062+2410 has increased during the period we have moni- 
tored it (see Figure[3}. 

A similar uncommon increase of the radio flux density was 
detected in CRL 618 from 1975 to 1991 dKnapp et alj|l995l) . If, 
as observed in CRL 618, we assume a linear increase, we can fit 
the three points in Figure [3] and extrapolate that the ionisation in 
IRAS 18062+2410 started around 1991, confirming the rapid evo- 
lution previously obse rved in this object I Parthasarath v et al.l2000l : 
lArkhipova et alfcOQll) . 

The spherical model used in modelling this object allowed us 
to match the 2003 data, using the 2005 1.4 GHz value as an upper 
limit. We determine that the shell in this target is optically thick be- 
yond 5 GHz. Simil ar behaviour has been observed in CRL 618 also 
(Knapp et al. 1995). Further observations, in both the radio and at 
other wavelengths, will help us to determine whether these objects 
share other common features. 

3.3 Radio morphology 

Figure Q] shows the structures observed in the resolved targets. All 
except IRAS 18062+2410 show bipolar radio morphology, with 
clearly resolved peaks. IRAS 17381-1616 shows a hint of bipo- 
larity, althou gh the angular resolution is not adequate to show its 
morphology. lAaquist & Kwold dl99lh also observed a sample of 
young PN with high resolution in the radio and found that most of 
their targets had a bipolar structure. Such radio morphologies, to 



4 SUMMARY 

We have presented radio continuum observations of a small sam- 
ple of pre-PN. The observations were carried out to inspect both 
the morphology and the physical conditions of the ionised com- 
ponent in these nebulae. Compact planetary nebulae are, in prin- 
ciple, only one possible explanation for radio-frequency contin- 
uum emission from evolved stars, the others being stellar photo - 
sphere, chromosphere, and circumstellar dust (Knap rTet al.lll995l) . 
Our observations rule out the possibility of a stellar wind or dust 
origin for the emission, since the obtained spectra are not compat- 
ible with such hypotheses. The dust would exhibit an approximate 
black body spectrum (S v oc v 2 in the radio range), and a stellar 
wind would have a 0.6 spec tral index over a wide range of frequen- 
cies dPanagia&Fellil 19751) . 

We have modelled the observed spectra assuming that the 
emission arises in spherical shells around the central stars, with 
the density decreasing as r~ 2 . The values of electron densities and 
ionised masses found for our targets confirm their nature as young 
PN, since they fall within the ranges expected for such objects. The 
data seem to point to a somewhat steeper than expected slope be- 
yond 5 GHz but, by spectral fitting, we calculated spectral indexes 
that match, to within errors, the theoretical value of -0.1. 

IRAS 18062+2410 appears to be optically thick up to 8.4 GHz 
and also shows flux density variations on a time scale of a few 
years. Such variations are rarely observed in the post-AGB/young- 
PN transition phase and indicate that this star belongs to a group of 
objects having such properties as those observed in CRL 618. IRAS 
17423-1755 also seems to be optically thick at high frequency, but 
the complete spectrum is needed to confirm this result. Quite inter- 
estingly, IRAS 20462+3416, whose overall centimetre spectrum is 
very flat, has a larger flux density at 22.4 GHz than at lower fre- 
quencies. We speculate that this may be due to a contribution from 
cold dust. 

The high-angular resolution observations that we have pre- 
sented show that the ionisation in these targets has already involved 
the walls of the cavity around the central star and is not limited 
to the tenuous gas within the cavity. All of the targets have been 
at least partly resolved, except IRAS 18062+2410. Almost all of 
the envelopes show two peaks of emission in a somewhat extended 
nebulosity that can be explained as an opacity effect, possibly due 
to a circumstellar torus, or the onset of jets. Unlike the other targets, 
the envelope in IRAS 22023+5249 shows more of an elongation 
along an axis than two bright peaks. Such a morphology would be 
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better explained by the action of jets rather than the presence of a 
torus. 
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